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Abstract 

We study a model of ^-component complex fermions with a kinetic term that is second order 
in derivatives. This symplectic fermion model has an Sp{2N) symmetry, which for any N contains 
an SO(3) subgroup that can be identified with rotational spin of spin-^ particles. Since the spin-^ 
representation is not promoted to a representation of the Lorentz group, the model is not fully 
Lorentz invariant, although it has a relativistic dispersion relation. The hamiltonian is pseudo- 
hermitian, H< = CHC, which implies it has a unitary time evolution. Renormalization-group 
analysis shows the model has a low-energy fixed point that is a fermionic version of the Wilson- 
Fisher fixed points. The critical exponents are computed to two-loop order. Possible applications 
to condensed matter physics in 3 space-time dimensions are discussed. 
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I. INTRODUCTION 



A basic result in the quantum field theory of fundamental particles in 4 space-time di- 
mensions is that requiring Lorentz invariance for spin-| particles necessarily leads to the 
free Dirac lagrangian. Here "spin-|" refers to the 3-dimensional rotational subgroup of the 
Lorentz group. The Lie algebra of the Lorentz group is SU(2) ® SU(2), where one linear 
combination of the two SU(2) symmetries is identified as angular momentum; thus spin 
representations are promoted to Lorentz representations in a straightforward manner. 

This paper in part deals with the following basic question. Suppose one wishes to describe 
the quantum field theory of spin-| particles in a physical context where the dispersion 
relation happens to be Lorentz invariant, but unlike in fundamental particle physics the 
full Lorentz invariance is not evidently required. In particular, we have in mind systems in 
condensed matter physics where the relativistic dispersion relation arises as a consequence 
of special properties of the system, and what plays the role of the speed of light is some 
material-dependent velocity. For instance, the effective mass of electronic quasi-particles 
may go zero because one is near a quantum critical point or more simply because of band 
structure as in 2-dimensional graphene , and massless particles, like photons, usually 
require a relativistic dispersion relation. Under these circumstances, the question is whether 
such a quantum field theory must necessarily be that of a Dirac fermion. In the case of 
graphene, the fermionic quasi-particles do turn out to be described by the massless Dirac 
equation. The reason for this is not an intrinsic Lorentz invariance, but rather that the 
continuum limit of a tight binding model on a hexagonal lattice gives a hamiltonian that is 
first order in derivatives, and near the Fermi points the particles are massless. 

In this paper we will study an alternative to the Dirac theory. The model is built out 
of an iV-component complex fermionic field with a non-Dirac, two-derivative action with 
a Lorentz-invariant dispersion relation. This model has a symplectic Sp(2N) symmetry. 
If this symmetry is viewed as an internal symmetry, then the fields are Lorentz scalars 
and the theory is Lorentz invariant but with the "wrong" statistics. However, since the 
Lie group Sp(2N) has an 5*0(3) subgroup, we can identify the latter with rotational spin. 
Therefore our model can naturally describe spin-| particles. The fermionic statistics is 
then in accordance with the spin-statistics connection, which requires spin-| particles to 
be described by fermionic fields. However since the rotational spin-| is not promoted to a 
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representation of the Lorentz group as in the Dirac theory, our model is strictly speaking 
not Lorentz invariant; hence the terminology "semi-Lorentz" invariant. 

The most serious potential problem of this theory concerns its unitarity, and this is 
addressed in the present paper. We show that the hamiltonian is pseudo-hermit ian, 

W = CHC, (1) 

where C is a unitary operator satisfying C 2 = 1. This generalization of hermiticity was 



considered long ago by Pauli 



3[], and more recently by Mostafazadeh 4] as a way of explaining 



;he real spectrum and addressing the unitarity issue in VT symmetric quantum mechanics 
^,0]. The important point is that by suitably defining a C-dependent inner product, pseudo- 
hermiticity of H is sufficient to ensure a unitary (i.e., norm-preserving) time evolution. This 
is explained in section II. One should also point out that taking a non-relativistic limit in 
the kinetic term one obtains a perfectly unitary second-quantized description of interacting 
fermions. 

The identification of the 50(3) Lie sub-algebra of Sp(2N) with rotational spin is described 
in section III. There we also study the discrete space-time symmetries of time-reversal and 
parity, and show how the spin generators transform properly under them. 

Another possible signature of non-unitarity comes from studying finite-size or finite- 
temperature effects. We show in section IV that whereas imposing periodic boundary con- 
ditions leads to a negative coefficient in the free energy, correctly imposing anti-periodic 
boundary conditions, as is normally appropriate for fermions, leads to a positive coefficient. 
(In two dimensions this coefficient is related to the Virasoro central charge.) 

Symplectic fermions are interesting for potential applications to critical phenomena. First 
of all, in D = 3, since the group of spacial rotations is simply U(l), there are less constraints 
coming from Lorentz invariance. More importantly, unlike Dirac fermions, four-fermion in- 
teractions in D = 3 drive the theory to some novel low-energy fixed points that are fermionic 
versions of the Wilson-Fisher p| fixed points. The reason is simple: in three dimensions a 
symplectic fermion x nas classical (or mean field) scaling dimension 1/2, whereas Dirac 
fermions if) have dimension 1; therefore x 4 is a dimension-2 operator and is relevant in 
D = 3, whereas ip 4 is irrelevant. This was the original motivation for the work [8j, where 



it was attempted to interpret these fixed points at iV = 2 as examples of deconfined quan- 
tum criticality Q|. Although it remains unclear whether our fermionic critical point can 



3 



correctly describe deconfined quantum criticality as defined in |9|, the resolution of the 
unitarity problem as presented in this paper is sufficient motivation to analyze the critical 
exponents further. In sections V and VI we extend the analysis of ja] to two-loop order. In 
particular, we show that some of the critical exponents can be obtained by analytically con- 



tinuing known results for the O(M) Wilson-Fisher fixed point to M = —2N [23(|. We also 
calculate the critical exponents for composite bilinear operators, which to our knowledge 
have not been studied for the 0(M) fixed points \2.\ . 

The correspondence of our model with O(M) for negative M is merely formal and not 
expected to be valid for all physical properties. First of all, the symmetries of the models 
are different. It should also be clear from the fact that in applications to condensed matter 
physics, our model has a Fermi surface, etc. Some concrete distinctions in the finite-size 
effects are made in section IV. 

In section VII we speculate on some possible applications to 2 + 1 dimensional quantum 
criticality in condensed matter physics. In the broadest terms, at the fixed point our model 
can be viewed as a quantum critical theory of spinons, where the symplectic fermions are 
fundamental spinon fields. For N = 2 components, we discuss possible applications to quan- 
tum anti-ferromagnetism, where the magnetic order parameter n is a composite operator 
in terms of spinons n = x^&X- This compositeness is the same as in deconfined quantum 
criticality {J, however our model is different since it has no U(l) gauge field. We show that 
two-point correlation exponents (rj) are rather large compared to the bosonic Wilson-Fisher 
fixed point, and this is mainly due to the compositeness of ft. By treating both cases, we 
show this is true irrespective of whether the particles are bosons or fermions. 

Section VIII contains a summary of our main findings and some conclusions. 

II. PSEUDO-HERMITICITY AND UNITARITY OF COMPLEX SCALAR 
FERMIONS 

Let x denote an iV-component complex field and consider the following action in D = d+1 
dimensional Minkowski space: 



S x = I efxdt 



d^xd^X -m XX~4n g ( X x) , (2) 
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where XX — Y,a=iX a X a '■> and d^d^ = d" 2 — S 2 .. If % is taken to be a Lorentz scalar, then 
the model is Lorentz invariant. The above action has an explicit internal U(N) symmetry. 
In the next section we show that, if x is a fermion, then there is actually a hidden Sp(2N) 
symmetry. 

We wish to quantize this model with x taken to be a fermionic (Grassman) field. The 
conjugate-momentum fields are 7r x = —dix and vr^ = d t X- They obey the canonical anti- 
commutation relations 

{T (x, t), d tX p (x', t)} = -{x°(x, t), ^(x', t)} = iF*6W(x - x') . (3) 

The hamiltonian for this system is 

H = J d d x [dtxdtX + d x xd x x + ™ 2 XX + 4n 2 g [xxf] • ( 4 ) 

Note that because of the fermion statistics the interaction term vanishes for N = 1. 

Consider first the free theory with g = 0. Suppressing the component indices, the fields 
have the following mode expansions consistent with the equations of motion: 



Ak-x 



where uu^ = Vk 2 + m 2 and k ■ x = uj^ t — k ■ x. The extra minus sign in the b_ term in 
X is chosen so that the anti-commutation relations ([3]) lead to the standard non-vanishing 
canonical relations 

{&k+, bl + } = {6 k _, bl_} = *«>(k - k') , (6) 

where it is understood that the field operators belong to the same field component. Of 
course, the extra sign would be unnecessary if x was a bosonic field. 

With the above definitions the fields have the required properties under causality (see for 
instance [ll(), namely x^iu)} — {x a { x )i X^iu)} — and 



{X a (x), X P (V)} = 5 a(S [A+{x -y)- A + (y - x)) , (7) 

where 
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Since A + (x) depends only on x 2 and is well-defined for space-like separation (x 2 < 0), 
it follows that x( x ) an d x{y) anti-commute for space-like separated points x and y. This 
implies that the hamiltonian densities TC(x) and TC(y) also commute at space-like separation. 
We will return to the spin-statistics connection below. 

In terms of the momentum-space modes, the free hamiltonian is 

H = J d d k tu k (fl + b k+ - &k_&t_) = J d d k cu k (bt + b k+ + &L&k-) + const. (9) 

We can discard the infinite constant in the above equation, which is equivalent to normal- 
ordering the hamiltonian. Let us define the vacuum |0) as the state being annihilated by b±. 
As a result, all states have positive energy. The one-particle states are doubly degenerate: 

Hbl ± \0) = H\k,±)=u k \k,±). (10) 

For N — 1 the theory is manifestly invariant under a U(l) symmetry. The corresponding 
conserved current satisfying <9 M J M = is J M = i[(d^x)X ~ X^/iX]- The conserved charge 
Q = J d d x Jt(x) can be expressed as 

Q = J d d k u k (b[ + b k+ + &,<-&£_) • (11) 

With these conventions, the operators 6+ and b_ have charge Q — 1, whereas b^_ and b + 
have charge Q = —1. 

The usual spin-statistics connection (see for instance [ll|) is based on causality as de- 
scribed above, along with the requirement that the hamiltonian must be constructed out 
of x an d its hermitian adjoint in order for it to be hermitian. The latter is violated here: 
because of the extra minus sign in the mode expansion of x f° r the fermionic case, one sees 
that unlike for the bosonic case, x is n °t the hermitian adjoint of x- Ld; us introduce a 
unitary operator C satisfying C^C = 1 and C = C\ which is defined by the properties 
Cb±C = ±b± and Cb±C = ±6j_. Then the relation between x an d X can t> e expressed as 

X = C X f C. (12) 

Since (xxY = CxxC, the hamiltonian satisfies the "intertwined" hermiticity condition 

H^ = CHC. (13) 

The above is also true in the interacting theory since [(xx) 2 ]^ = C{xx) 2 C- (Note that for the 
free theory in momentum space, C actually commutes with H as can be seen from eq. ©, 



and hence the above equation is trivially satisfied; however, this is no longer the case in the 
interacting theory.) 

n 

Hamiltonians satisfying eq. (TTHT) were considered long ago by Pauli |3[ and more recently 
by Mostafazadeh |4| in connection with VT symmetric quantum mechanics We 
will follow the previously introduced terminology and refer to such hamiltonians as C- 
pseudo-hermitian. Quantum mechanics based on pseudo-hermitian operators has some very 
desirable properties, which parallel the standard ones. First of all, as we now explain, a 
pseudo-hermitian hamiltonian can still define a unitary quantum mechanics if one defines 
the inner product appropriately. Specifically, consider a new inner product defined as 

ww) c = mew ■ (14) 

Then probability is conserved with respect to this modified inner product, i.e., the norms of 
states are preserved under time evolution: 

W(t)\il)(t)) c = W\e iHH Ce- lHt \^} = {i)'\Ce iHt C 2 e- lHt \i)) = (^) c . (15) 

Pseudo-hermiticity also ensures that the eigenvalues of H are real. To see this, let \i[>e) 
denote an eigenstate of H with eigenvalue E. Then 

(E - E*) (V*>c = (iPe\{CH - H^C)\^ E ) = . (16) 

Therefore eigenstates of H with non-zero C-norm necessarily have real energies. Other 
properties are proven in 4|. 

In the sequel it will be convenient to define the pseudo-hermitian adjoint A' c of any 
operator A as the proper hermitian adjoint with respect to the C-inner product: 

mA\^r c = (i/>\Ai'W)c, (17) 

which implies 

A^ = CA ] C. (18) 

The pseudo-hermiticity condition on the hamiltonian then simply reads = H. One can 
easily establish that this pseudo-hermitian adjoint satisfies the usual rules, e.g. 

(AB) tc = B^A U , 

(aA + bB) u = a*A^ + b*B u , (19) 

where A, B are operators and a, b are complex numbers. 



III. Sp{2N) SYMMETRY AND ROTATIONAL SPIN 



The action (j2j) has an explicit U(N) symmetry irrespective of whether x is bosonic or 
fermionic. If \ is bosonic, then the action expressed in terms of real fields has an 0(2N) 
symmetry. On the other hand, if x is a fermionic field, then the symmetry is Sp(2N). Even 
for N = 1, the Sp(2) symmetry of the fermionic theory is larger than the U(l) symmetry of 
the bosonic theory. To manifest this symmetry explicitly, let us express each component \ a 
as 

X a = ^ (V? + irg) , T = ^= K - «£) • (20) 
We now introduce the 2x2 anti-symmetric matrix e = ( ^ ) and the 2N x 2N matrix 
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= e (g) Itv- Arranging the real fields rj" into a 2N vector r\ = (ri\,ril,rif,ri%, . . . ,77^,77^), 
one has 

XX = 2 V T £nV ■ ( 21 ) 

This bilinear form has the symmetry rj — > Mr], where M ' < \ M = €n- This is the defining 
relation for M to be an element of the group Sp(2N). 

We will need the Lie algebra of Sp(2N). Let M = e x , in which case the above relation 
implies X T e^ = — e^X. A linearly independent basis for X is then {1 <g> A, o x £g> t x ,a y ® 
t y , <J z ®t z }, where A is an iV x anti-symmetric matrix, Oi are the Pauli matrices, and ti are 
N x N symmetric matrices [3]. For any iV the algebra Sp(2N) has an 5*0(3) sub-algebra 
generated by o"j <g> 1^, which can in principle be identified with spin. It also has an SO(N) 
sub-algebra generated by the matrices A, and an SU(N) sub-algebra generated by I® A and 
o z ® t z , where t z is traceless. The Lie algebra Sp{2) is equivalent to SO(3) = SU(2). Note 
that the N = 2 case therefore has two different 50(3) sub-algebras that could potentially 
be identified with spin. 

A. Canonical quantization and pseudo-hermiticity 

For simplicity, let us specialize to the free theory with N — 1 component (the generaliza- 
tion to iV 7^ 1 is straightforward). Then the action takes the form 

S = Jd d x dt (d^id^ - m 2 ^) • (22) 



S 



The canonical momenta are tti = dtt]2 and 7C2 = —dtTji, which leads to the equal-time 
anti-commutation relations 

Mx^^x',*)} = -{^(x,*),^!^',*)} = ^(x-x'). (23) 

The pseudo-hermiticity of the hamiltonian exactly parallels the discussion in section II. 
If one expands the fields as 

* w = /(2^7sK- e "" tI + a -^)- (24) 

then eq. (1231) implies 

{a^,al_} = {a k+ ,4'+} = ^ d) (k- k') (25) 
In terms of these modes, the hamiltonian is 

H = J d d ku* (a k+ a k+ - a k _4_) • (26) 

The relation between 7] 2 and 771 is r\% = Crj{C, where C flips the sign of a_ and aL, i.e. 
Ca±C = ±a± and Ca±C = ±a±. The hamiltonian is pseudo-hermitian, = CHC, in 
both the free and interacting theory. Thus, for the reasons explained in section II, it defines 
a unitary time evolution. 



B. Sp(2N) conserved charges 

For N = 1 the Sp(2) symmetry is 

V =^^e x V , (27) 

where X = a ■ a with arbitrary parameters c^. The conserved currents following from 
Noether's construction read 

J l = ~2 (Vid^ + md^Vi) , Jt = -iVid^Vi , J? = ^2^2 • (28) 

Using the equations of motion and fermion statistics (rjf = 0) one readily verifies that 
<9 M J M = 0. The conserved charges are then defined as usual as S = J d d x J t . In terms of the 
creation and annihilation operators, they take the form 

S z =^ J d d k (a k+ a k+ - a k _a k -) , S + = J A4 + a k _ , S~ = J Aa^a^ . (29) 



As expected, these charges satisfy the Sp(2) = £0(3) algebra [S^jS^] = ±S^, [S + ,S~] = 
2S Z with the identification S* 1 * 1 = S x ± iS y . 

Using the pseudo-hermiticity of the hamiltonian, one finds that the above conserved 
charges have the pseudo-hermiticity properties (S z )' c = S z , (5 ,± )^ c = —S T . This implies 
that pseudo-hermitian conjugation of the Sp(2) generators is an inner-automorphism of the 
algebra: 

(S) u = e i7r& Se-™ s * . (30) 

C. Identifying the spin 

The spin-statistics connection requires that particles with half-integer spin under rota- 
tions be quantized as fermions. Since the Lie algebra Sp(2N) has an 50(3) sub-algebra 
generated by a® In, h is natural to try and identify this 50(3) sub-algebra with spin (i.e., 
spacial rotations). 

Again for simplicity let us consider N = 1 component symplectic fermions. The one- 
particle states |k, ±) = a k± |0) of energy u; k have spin S z |k, ±) = ±| |k, ±), so that a\_ and 
al create spin-up and spin-down particles, respectively. A further check of this identification 
comes from considerations of time-reversal symmetry, to which we now turn. 

D. Time reversal and parity 

Let T denote the time-reversal operator. Since T is anti-linear, it can be written as 
T = UK, where U is unitary and K complex conjugates c-numbers: Kz = z*K. Consider 
spin- 1 particles, where spin is represented by the Pauli matrices a. Since spin is odd under 
time reversal, T<rT _1 = —a, which implies Ua*U* = —a. The well-known solution to this 
equation is U — o y . Since time reversal also flips the sign of momentum, we are led to define 

Ta^T' 1 = ±i a_ kjT , Ta k ± T _1 = =p al k>T . (31) 

As is well known, due to the anti-unitarity, T 2 = — 1 on one-particle states of spin ~. 

Using the above transformations in eq. ( |29i) . one sees that the Sp(2) generators have the 
correct transformation properties to be identified as rotational spin: 

TST- 1 = -S. (32) 
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From the form of eq. (126]) it follows that the hamiltonian is invariant under time- reversal, 
i.e. THT- 1 = H. 

On the modes, parity simply flips the sign of momenta, i.e. Va^V = a_k,± ; an d similarly 
for a±. The hamiltonian is thus also invariant under parity. 



IV. FREE ENERGY AND FINITE SIZE EFFECTS 



Finite-size effects are another probe of the unitarity of a theory. Let us therefore consider 
our model embedded in the geometry of ci-dimensional flat space with periodic time described 
by a circle of circumference /3, i.e., i? d ®S' 1 . We will use the language of quantum statistical 
mechanics and identify f3 = 1/T with T being the temperature. The <i-dimensional volume 
of R d will be denoted as V. 

In Euclidean space the action is 



J d D x [dxdx + m 2 xx + ^ 2 g {xxf\ , (33) 



where D = d + 1 is the Euclidean space-time dimension, d 2 = J2iLi d 2 , and as before x is 
an iV-component complex fermion field. In order to make certain arguments in the sequel, 
let us introduce an auxiliary field u(x) and consider the action 

,2" 



S x , u — J ' 



d D x 



dxdx + (m 2 + 2ixyfg u) XX ~ "j 



(34) 



from which the original action S x is recovered when the field u is eliminated using its 
equations of motion. Since x now appears quadratically, one can perform the functional 
integral over it to obtain 

Z = J VxDxVue~ Sx ' u = J D«e- Scff , (35) 

where 

/u 2 
d D x — , (36) 

and we have used the identity logdetA = TrlogA Note that if x were taken to be a 
complex bosonic field, then the functional integral would give 1/ det A rather than det A, 
which amounts to the replacement N — > —N in Seff. This suggests that some physical 
quantities in the symplectic fermion model can be obtaining by flipping the sign of N in 
its bosonic counterpart. However, we now demonstrate that this not correct for all physical 
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quantities; in particular, such a replacement does not hold for the free energy when one 
takes into account the proper boundary conditions. 

For the remainder of this section we will consider the non- interacting theory (g = 0). In 
order to be able to perform the integrals and to compare with some known results, we also 
set the mass m to zero. In the free theory the field u decouples, and the functional integral 
over u just changes the overall normalization of the partition function Z. Discarding this 
overall factor one obtains Z = exp [N Tr log(— d 2 )]. The free energy density T = —T log Z/V 
is then simply T = — ^y- Tr log(— d 2 ). With the Euclidean time compactified on a circle of 
circumference /3, the time component of the momentum is quantized, ko = (2iru//3), where 
v is a Matsubara frequency. The functional trace is then 

Tr log(-d 2 ) ^ E /(g,log[k 2 + (2nv/{3) 2 } . (37) 

In the Euclidean functional integral approach to finite temperature, one is required to 
impose periodic boundary conditions for bosons and anti-periodic boundary conditions for 
fermions. In order to illustrate an important point, let us first consider \ to be periodic, i.e. 
v an integer. It is a well-known identity that 

£ log [k 2 + (27rz///3) 2 ] = (3u k + 2 log(l - e~^) , (38) 

where = Vk 2 for m — 0. The first term above gives a temperature-independent contri- 
bution to the free energy, which we can discard by defining T such that it vanishes at T = 0. 
The result is 

^=-2JVT|^log(l- C -^). (39) 
In analogy with black-body formulas in four dimensions, let us define a coefficient Co through 

^-,»T» (40) 



7T 



where ( is Riemann's zeta function. The above normalization is such that cp = 1 for a 
single free massless boson. Performing the integral in eq. ( |39l) one obtains 



c D = -2N (periodic b.c.) (41) 

in any dimension D. The negative value of cd is normally a sign of non-unitarity. In two 
dimensions, for unitary theories with zero ground-state energy, C2 is the Virasoro central 
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charge Cv; r of the conformal field theory 13J, and Cyir 2 IS the usual result for a single 



symplectic fermion [14j . (For a precise, general relation between C2 and c v ; r see the end of 
this section.) Note also that cd = — 2N is simply the N — > —N result for iV free complex 
massless bosons. This is to be expected, since we have computed it using the periodic 
boundary conditions appropriate to bosons. 

The result ( l4Tj) is incompatible with the spectrum of particles computed in section II. 
Clearly this is due to having taken the wrong boundary conditions for the fields. For anti- 
periodic boundary conditions, v is half-integer, and one has 

log [k 2 + (2nu/(3) 2 ] = (3u k + 2 log(l + e"^) (42) 

vez+i/2 



F = -2NT I -?— d log(l + e'^) . (43) 



instead of fl3"g|) . leading to 

d d k 

(2tt) c 

It is clear from the above expression and basic results in quantum statistical mechanics 
that this result corresponds to 2N free fermionic particles with one-particle energies uj^, 
consistent with the quantization in section II. Performing the integral one finds 

c D = 2N ^1 - ^75^1^) (anti-periodic b.c.) (44) 

in D dimensions. The thermal central charge cd is now positive and consistent with a 
unitary theory; in fact, it is the same as for 2N real Dirac fermions. 

A few additional remarks clarifying the well-studied 2D case are in order. Consider the 
first-order action 

S b/C = J d 2 x(bd-c + bd z c) } (45) 

where z and ~z are Euclidean light-cone coordinates. Let us define Lorentz "spin" with 
respect to Euclidean rotations, such that for a holomorphic field ip s (z) of spin s 

ij s (e 2irr z)=e 2 ™ s ^ s (z). (46) 

With this convention, usual Dirac fermions have spin s = ±|. Note that parameterizing 
s = i9/ (27r) implies that the spin s is defined modulo = 2tt. Let us assign the following 
spins to the b and c fields: 

spin(6, c) = + ^, -s + ~) , spin(6, c) = (s - i, -s - ^) . (47) 
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Then the Virasoro central charge is known to be C2 = 1 — 12s [15(. Identifying 

b = d z \ ] , ckc = ckx , b = ckx ] , d z c = d z \ , (48) 

then the above first-order action is equivalent to our symplectic fermion action with x ~~ * X ■ 
The above identifications are consistent with spin(x) = — spin(%T) = | — s (with defined 
modulo 27r as above). The usual correspondence between symplectic fermions and first-order 
actions is based on letting Xi X* have spin 0, which implies s = | and C2 = —2 [l4]. However 
another choice is spin(x) = — spin(x^) = ~, which gives s = and C2 = 1, as we found above 
in the 2D case. 

Another check in two dimensions goes as follows. The thermal central charge C2 in this 
section is known to be related to the Virasoro central charge c v j r by the formula c e g = c v ; r — 
24A min , where A min is the minimal conformal scaling dimension. In the twisted (Ramond) 



sector of the symplectic fermion, t 
with dimension A m i n = —1/8 



die ground state is known to correspond to the twist field 
Since twist fields modify boundary conditions from 
periodic to anti-periodic, a consistency check is that the value C2 = 1 in eq. (jUJ) for N — 1 
should correspond to c e g with c vir = —2 and A min = —1/8, and indeed it does. 



V. RENORMALIZATION GROUP AND CRITICAL EXPONENTS 

We study the interacting critical point of the symplectic fermion theory described by 
the Euclidean action in eq. (|33|) . using a position-space approach based on the operator 
product expansion (OPE). In the following section we will present an alternative derivation 
of the critical exponents (extended to two-loop order) using a technique based on Feynman 
diagrams familiar from quantum field theories for elementary particles. 

Consider a general Euclidean action of the form 

S = S + 4tt 2 2 J d D x 0{x) , (49) 

where So is conformally invariant, g is a coupling, and O a perturbing operator. For our 
model, So is the massless free action and O = (xx) 2 - To streamline the discussion, let [X] 
denote the scaling dimension of X in energy units, including the non-anomalous classical 
contribution which depends on D. An action S necessarily has {SJ = 0. Using [<i D x] = —D, 
the classical dimensions of the fundamental couplings and fields are determined to be [x] = 
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(D — 2)/2, [m] = 1, and [g] = 4 — D = e. Let us therefore define the quantum corrections 



to the scaling dimensions 7 X and •y m as [25j 



ixl = ^j^ + l x , [mj = l + j m . (50) 
At the critical point, r y x determines the two-point function of the x fields via 

(X f (x)x(0))~ |x p_ 2+27x - (51) 

The anomalous dimension 7 m can be used to define a correlation-length exponent v. At the 
critical point, the correlation length diverges as m — > 0, i.e. £ ~ mT 2u . Using the fact that 
[£] = — 1, one has — 2u = [£]/[m], which implies 

1 



" 2(1 +7m)'" (52) 

The lowest-order contributions to the /3-function and the critical exponents are easily 
calculated in position space. At first order in the ^-expansion, the OPE coefficients can be 
computed in four dimensions. Consider first the /3-function. Since O is a marginal operator 
in D = 4 (classically [O] = 4), the OPE gives 

0{x) 0{y) = ° 0{y) + ... (53) 

47T 4 |X — |/| 4 

for some coefficient C. Consider now the correlation function (X) for arbitrary X to second 
order in g: 

(X) = (X) - 4n 2 g J d'x (X O(x)) + \ {^ 2 gf J d'x J d*y (X 0(x) O(y)) + ..., (54) 

where the subscript indicates that the correlation function is computed with respect to 
the free action So. Using the OPE (153]) in the above expression along with J a d A x/x A = 
— 27r 2 loga, where a is an ultraviolet cut-off, one finds 

(X) = (X) - A7i 2 {g + Cg 2 log a) J d 4 x (X 0{x)) + ... . (55) 

The ultraviolet divergence is removed by letting g — ► g(a) — g — Cg 2 log a. This leads to 

P(g) = - 7 ^ = -eg + Cg 2 + ... J (56) 
a log a 

where the leading term comes from the classical dimension of g. (Our convention for the 
sign of the beta-function is as in high-energy physics, where increasing a corresponds to a 
flow toward low energy.) 
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The above calculation easily generalizes to actions of the form 

3 = S + 4vr 2 9a [ d D x O a {x) , (57) 

which typically arise in anisotropic versions of our model. If the perturbing operators satisfy 
the OPE 

O a (x) O b (y) = £ An , lx c _ yl , OM + ..., (58) 
then the corresponding /5-functions are 

Pa(g) = -Sg a + Y, C a9o9c + -- - ■ (59) 

b,c 

Let us return now to our model. Using the OPE results 

Xi(x)xM = -Xi(x)Xj(0) ~ -4^^ (60) 

valid in four dimensions to evaluate (1551) for O = (xx) 2 , one finds C = A — N which leads to 

P(g) = -sg + (A-N)g 2 + ... . (61) 

The model thus has a low-energy fixed point at ~ ej (4 — N). 

Consider next the anomalous dimension of the symplectic fermion fields and of composite 
operators built out of these fields. Let $(x) denote a field having the following OPE with 
the perturbation: 

0{y) *{x) = y[4 $(y) + ■ ■ ■ (62) 

for some coefficient B. Then to first order in g 

($(0)) = ($(0)) - 4n 2 g J d 4 y (0(y) $(0)) + . . . 

= (1+Bg log a) ($(0)) + . . . « a B ° ($(0)) . (63) 

This implies an anomalous contribution 7$ to [<&] given by 7$ = + . . .. For the operator 
$ = XX-, the OPE result (l60l) implies B = 1 — N. Finally, using 2[m] + = D, one has 



1 ^-1 

lm = ~-lxx = ^—9 + ■■■ ■ (64) 
In the sequel it will also be of interest to consider other fermion bilinears of the form 

n T (x) = xrx, (65) 
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where r is a traceless matrix. This operator does not mix with XX i n the OPE f|62|) and has 
an independent anomalous dimension 7„. Repeating the above computation, one finds that 
because of the tracelessness of r there is no contribution to this order proportional to N, 
i.e., B — 1, and this leads to 

7„ = g + . . . . (66) 

To the order we have computed so far, our results for the /3-function and the anomalous 
dimension 7 m are the same as for the O(M) Wilson-Fisher fixed point with the substitution 
M = —2N. This evidently follows from the auxiliary-field construction in section IV: bosons 
versus fermions differ by the overall sign of the logarithm of the determinant, which amounts 
to N — > —N in the effective action S e g in eq. fl36l) . The factor of 2 in M = —2N comes 
from the fact that x is a complex field, whereas the O(M) model is formulated in terms 
of M real fields. Though one may worry that this equivalence with exponents of 0(—2N) 
may be spoiled at higher orders for certain operators whose correlation functions cannot be 
computed from S e g, we verify in the next section that the equivalence persists to two-loop 



order. We thus expect it to hold to all orders in perturbation theory 



26|]. 



VI. TWO-LOOP RESULTS 

The simple position-space method of the last section does not extend straightforwardly 
to higher orders. Here we describe an alternative calculation of the /3-function and the 
anomalous dimensions using Feynman graphs. 

We consider the action (T5]) of an iV-component symplectic fermion x in Minkowski space 
and express it in terms of bare parameters mo and go and unrenormalized fields Xo '■ 



Sy^ — Id X 



duXo d "Xo - ml XoXo ~ ^ 2 9o(x Q Xo) 2 ■ (67) 



The momentum-space Feynman rule for the four- fermion vertex with incoming fermions x", 
X 13 and outgoing fermions x 13 is ( — ^ 2 Wo) if ot ^ /3, while it vanishes for a = /3 due to the 
anti-commuting nature of the fields. The momentum-space propagator for the fermion x a is 
diagonal in component indices and given by the ordinary Feynman propagator for a scalar 
field, i/(p 2 — tuq + iO). The mass term is kept in our calculations as an infrared regulator. 

The mass dimensions of the field and coupling are [xo] = (D — 2)/2 and {g } = 4 — D. 
We work in dimensional regularization and define the renormalized coupling g(fi) through 
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FIG. 1: Four-fermion vertex function relevant to the calculation of the (3- function. 



go = ~p 4 ~ D Z g ([i) g(fi), where /I 2 = /i 2 e 7B /(47r). The scale \x acts as an ultraviolet regulator in 
momentum space. Our renormalization factors will be defined using the modified minimal 
subtraction (MS) scheme in D = 4 — e space-time dimensions. 



A. Calculation of the /5-function 

Standard field-theory arguments can be used to show that the /3-function of our model is 
given by 

i 3(g,D) = ^- = i3(g) + (D-A)g, (68) 
ain/i 

where 

(69) 

is independent of D. The quantity Z^ denotes the coefficient of the 1/e pole in the Laurent 
expansion of the renormalization factor near e — 0. Throughout, we denote g = g(/i) unless 
otherwise noted. 

The /^-function is obtained from the four-fermion vertex function with ingoing and outgo- 
ing component indices a (3, shown in Figure [TJ For simplicity we set the external momenta 
to zero. The tree-level contribution to the vertex function is given by the elementary vertex 
shown in Figure EJ At one-loop order there is a single loop topology but three different 
contractions of indices, depicted in Figure [31 which yield a multiplicity factor of (4 — N). 
At two-loop order there are three loop topologies, whose multiplicities can be obtained by 
analyzing the various possible contractions. The relevant diagrams and their group-theory 
factors are depicted in Figure HI The two-loop scalar integrals required for our calculation 
can be obtained from 171 |. 



Adding up the results for the various diagrams we obtain the bare vertex function. We 
then multiply this result by Z 2 to account for wave-function renormalization, and substitute 
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FIG. 2: Tree- level contribution to the vertex function. 






=> 



4-N 



FIG. 3: One-loop contributions to the vertex function. The three graphs give rise to the same loop 
topology, shown on the right. The group-theory factors of the individual diagrams are — (N — 2), 
1,1, where the minus sign of the first graph results from the closed fermion loop. 



Z m 2 m' 2 and go = /i 4 D Z g g to implement mass and coupling-constant renormalization. 



The renormalization factors Z x and Z m 2 are determined from the calculation of the fermion 
self-energy in the next subsection. By requiring that the renormalized vertex function be 
finite, we extract 

"(4 -iV) 2 3(3iV-7)" 



1 



N 



■9 + 



+ 



4s 



9 2 + 0(g* 



(70) 



From eq. ( 1691) we then obtain for the /3-function 



m = (4 - N) g> + 3(3A ! 7) g' + 0{g*) 



(71) 



The D-dimensional /3-function in ( 168|) has a non-trivial fixed point at positive coupling 
given by the solution to the equation f3(g*)/g* = 4 — D = e. At second order in the 
e-expansion, we find 



9* 



3(7 -3A0 2 3 , 

+ 777-, 77777 £ +0 (e 3 ) 



N 2(4 -N) 3 



(72) 
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FIG. 4: Two-loop topologies for the vertex function. Each topology receives contributions from 
several diagrams. The combined group-theory factors are listed below each graph. 

B. Calculation of the self-energy 

Next we need the anomalous dimension of the fermion mass and field. They follow from a 
two-loop calculation of the self-energy £(p 2 , m 2 ,) in the vicinity of the mass shell (p 2 = m 2 ). 
Here m and m are the bare and renormalized mass parameters, respectively. The relevant 
relations are 

m 2 = m 2 + £(m 2 , m 2 ) , Z~ x = 1 - S'(m 2 , m 2 ) , (73) 

where the prime denotes a derivative with respect to the first argument, p 2 . Inserting here 
777.Q = Z m 2m 2 , one finds for the renormalization factors 

Z m * = 1 - ^'Zf™ 2 ) ; Z" 1 = 1 - £'(m 2 ,Z m2 m 2 ) . (74) 

At one-loop order there is a single tadpole graph to evaluate, while at two-loop order we 
have a double tadpole diagram and the sunrise diagram, see Figure El After accounting for 
coupling-constant renormalization using eq. (I7D1) . we obtain 



1-N 



3^ 



Z x = l-^ r g 2 + 0(g i ). (75) 
The anomalous dimensions of the mass and field are given by 

7a = _„^2 1= _i d ^l_ (76) 

7m 9 dg ' Jx 2 dg • [ ' 



where in the second relation we take into account that xo — y Z x X- We find 

5 



-y m > = 2'Ym=(N-l)g(l--g > )+ 0(g 3 ) , 

l x = il ^ 1 9 2 + 0(g'). (77) 
20 
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FIG. 5: One- and two-loop diagrams contributing to the fermion self-energy. 



The anomalous dimension of the field starts at two-loop order. Instead of the anomalous 
dimension of m 2 one could compute the anomalous dimension of the fermion bilinear xXi 
which is given by j% x = -^ m %. 

Evaluating our expressions at the fixed-point value of the coupling yields 



7x 



(N 



2(4 - N) 
(1-N) 



22- 13iV 
+ 4(4 - iV) 2 

e 2 + 0(e 3 ). 



(78) 



3(4 — N) 2 

As a crosscheck of our results, we note that the two-loop expressions for the /3-function 
and anomalous dimensions obtained in eqs. fTTTl) and fl77j) go over to the corresponding 
results of O(M) scalar field theory (see e.g. [18|]) if we identify M = —2N. Likewise, the 
fixed-point of the fermionic /3-function in eq. f!72|) is related to the Wilson-Fisher fixed point 
[3] by the same replacement rule. In some sense, our symplectic fermion theory may thus be 
considered as an analytic continuation of scalar O(M) field theory to negative M. However, 
as emphasized in section IV, this simple correspondence is not expected to hold for all 
physical quantities. 



C. Renormalization of the composite bilinear n T 

The diagrams contributing to the renormalization of the composite operator n T = x T X 
defined in eq. (1651) can be obtained by inserting this operator into the one- and two-loop 
graphs for the fermion self-energy shown in Figure [5l In the evaluation of these graphs it is 
important to use that the matrix r is traceless. Writing the bare current as (n T ) = Z n n T , 
we obtain 

q f 3 - N 5 - N\ „ o N 
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The anomalous dimension of the current is thus 



^ — N 

in = g-- 1 -g 2 + 0(g 3 ). (80) 



At the fixed point, this yields 



e (2-A0(ll + AQ ^ 2 
ln = A^N + 4(4-iV)3 £ + ° {£) - (81) 

This result will become important for the discussion in the following section. 



VII. POSSIBLE APPLICATIONS 

In this section, we speculate on possible applications of the above results. The most 
interesting context is quantum criticality in d = 2 spacial dimensions (i.e., D = 3). In the 
broadest terms, since the particles have spin ~, our model can describe a quantum critical 
theory of spinons. 

Whereas the Mermin- Wagner theorem rules out continuous phase transitions at finite 
temperature in d = 2, zero temperature quantum phase transitions continue to be of great 
interest. The best studied example is the quantum phase transition in 2d Heisenberg mag- 



nets, which is in the universality class of the 0(3) Wilson-Fisher fixed point [19|, |20| . One fea- 
ture of our model is that it can describe quantum critical points wherein the magnetic order 
parameter n is a composite operator in terms of the more fundamental fermion fields. This 
could in principle have applications to quantum phase transitions in the anti-ferromagnetic 
phase of Hubbard-like models, where the magnetic order parameter is a bilinear in the elec- 
tron fields. As stated in the introduction, if such electrons were described by the Dirac 
theory, the four-fermion interactions are irrelevant and do not generally lead to a low-energy 
interacting fixed point. That this is different in the case of symplectic fermions was the 
primary motivation for our work. 

Let us first review the definitions of the exponents for the usual Wilson-Fisher fixed point. 
The order parameter is an M-component real vector n with action 



S- 



WF 



J d D x dn ■ dn + m 2 n ■ n + A (n ■ n) 2 . (82) 



Some of the exponents are defined with respect to perturbations away from the critical point. 
Namely, consider 

S = S* + J d D x (tO £ + B-n), (83) 
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where S* is the critical theory, O e is the "energy operator", and B a magnetic field. For 
applications to classical statistical mechanics in 3D one identifies O e = n-n and t = T — T c , 
so that t oc m 2 . The usual definition of the correlation- length exponent u e via £ ~ t~ Ve then 
leads to 

1 _ 1 

The second fundamental exponent is related to the scaling dimension of n. It is conven- 
tional to parameterize this with 77 in the form 

D- 2 77 



n 



+ 77 , (85) 



2 2 

so that the two-point function at the critical point scales as 

*(°)> ~ j^j]^ ■ (86) 

The convention for the leading contribution to \n\ comes from the action (182]) . which implies 
that classically n has dimension (D — 2)/2. The parameter 77 is then given as 77 = 27^, where 
7^ is the quantum anomalous correction to the scaling dimension of n. 
The third exponent characterize the one-point function of n via 

(n) ~ iP ~ (87) 

which leads to 

"=B-t( d - 2+ '>- (88) 

Treating i? as a coupling gives [5] + [n] = -D, from which it follows that 

m D+ 2-„ 

[n] + ,SJ ' 

For the remainder of this section we consider the special case where N = 2. In discussions 
of deconfined quantum criticality 9], for the case of 0(3) symmetry, the 3- vector n is 
represented as a bilinear in "spinon fields" Xi n = X where here x is an = 2 
component complex field and a are the Pauli matrices. Note that n is an example of 
an operator n T defined in (165j) . For N = 2, the above representation of n is consistent 
irrespective of whether x is a boson or fermion, so we will treat both cases in parallel. If n is 
a rotational 3- vector, then \ is a spin-| doublet. This identification of spin is different than 
in section III, and is instead based on the SU (2) sub-algebra generated by {1 <g> A, a z ® t z }, 
where t z is traceless and symmetric (see section III for notations). 
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Simple considerations based on the renormalization group point to a possibly special role 
played by the N = 2 fermion theory. Suppose a model formulated in terms of an 0(3) n field 
is asymptotically free in the ultraviolet region. Then the coefficient of the free energy cp 
described in section IV equals 3 for any D. On the other hand, free x fields give C3 = 3N/2 
by formula (j44j) in 3D. which is the same as for the 0(3) n theory when N = 2. Therefore, 
the free energies match up properly in the ultraviolet for an 0(3) n field and an N = 2 
symplectic fermion. 

The model proposed in Q| is based on the CP 1 representation of the non-linear 0(3) 
sigma model, which involves an auxiliary U(l) gauge field. The CP 1 model is then modified 
by relaxing the non-linear constraint n ■ n = x^X = 1 an d making the gauge field dynamical 
by adding a Maxwell term, effectively turning the model into an abelian Higgs model. One 
appealing feature of this model is that because of the equivalence of CP 1 and 0(3) non- 
linear sigma models (at least classically), without the Maxwell term one has an explicit 
map between the non-linear n field and x-field actions. Though this model is a natural 
candidate for a deconfined quantum critical point, unfortunately the fixed point is difficult 
to study perturbatively, so it has not been possible to accurately compare exponents with 



the simulations reported in 



21 



22|]. 



Let us broaden the notion of deconfined quantum criticality to simply refer to a quantum 
critical point for an 0(3) vector order-parameter n, where n is composite in terms of spinon 
fields x- The fields x are interpreted as the fundamental underlying degrees of freedom, and 
the critical theory S* in eq. (1531) is the critical theory for x- The critical exponents v, 77, (3, 
and S defined above are then related to scaling dimensions of composite operators in the 5* 
theory. 

The anomalous dimension of n in the epsilon expansion follows from the results in sec- 
tion VI and will be denoted 7^ in what follows. Let us first consider the case where x is a 
bosonic field. As explained above, the exponents for bosonic verses fermionic theories are 
simply related by N — > —N. Specializing eq. (181j) to N — —2 in 3D (e = 1) one finds 
7^ ps 0.21 for bosonic x- This is quite large compared to the analogous result 7^ ps 0.02 for 
the 0(3) Wilson-Fisher fixed point. (The latter is well-known; as explained in section V, it 
corresponds to 7 X at N = —3/2.) For x a fermion, interestingly the 0(e 2 ) correction to the 
leading one-loop result vanishes for iV = 2, so that 7^ = 1/2. We can give an alternative 
estimate of 7^ by simply substituting g* = 11/16 from eq. (1721) into (IHOl) without expanding 
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in e. This gives 7^ ~ 0.33 for fermionic x- F° r both the bosonic and fermionic cases, the 
largeness of 7^ compared with the usual 0(3) fixed point is due to the compositeness of 
n. If we naturally identify O e with XX-, then the correlation-length exponent u E is given by 
eqs. ( 1521 and ( I7H1) evaluated with N = ±2. This leads to 

v e « 0.75 (bosonic) , v e « 0.42 (fermionic) . (90) 

There are at least two difficulties encountered if we attempt to compare with existing 



numerical simulations, such as those in 



21 



221 ]. The main one is that the simulations 



are performed with an action for the n field or for local lattice spin variables Si with a 
Heisenberg-like hamiltonian, rather than with fundamental spinon degrees of freedom x, 
and we do not have a direct map between the two descriptions. In particular, in a theory 
with fundamental x fields and the compositeness relation n = x <?Xi since x has classical 
dimension (D — 2)/2, one has 

[n] = D - 2 + ln . (91) 

Comparing with eq. fl85|) one finds rj = D — 2 + 27^ rather than the usual relation rj = 27^. 
In other words, ri now contains a purely classical contribution of D — 2, and this was used 

Q 

to argue that the r\ exponent was large in [9J. On the other hand, simulations based on 
n- field actions effectively force the classical contribution to \n\ to be (D — 2)/2 as in the 
Wilson-Fisher theory, suggesting that simulations measure 77 = 27^. In the fermionic theory, 
support for this idea comes from the fact that the two lowest orders of the e expansion give 
7^ = 1/2 in D = 3. 

The other difficulty is that, unlike for temperature phase transitions where t = T — T c , 
in the context of zero temperature quantum critical points it is not obvious what plays 
the role of the parameter t, or equivalently, the energy operator O e that determines the 
correlation- length exponent. In this context, t is a parameter in the hamiltonian that is 
tuned to the critical point. The most natural choice is O e = XX which implies t = m 2 and 
leads to eq. fl90l . However, another possibility could be t = m, which would lead to twice 
the values in eq. (1911 . corresponding to v e = 1/(1 + j m ). 

The above difficulties prevent us from establishing a definite connection with the simu- 
lations in 2l|, |22). In fact, the exponents for deconfined quantum criticality are currently 
controversial, since the two above works disagree strongly on the value of the exponent 77. 
However, we point out that if one identifies 77 = 27^, then our computed exponents are not 
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inconsistent with some of the exponents in 
sive. More specifically, the work 



2ll . |22j , although the comparison is not conclu- 
2l| reports r] « 0.6-0.7 and v = 0.8-1.0. On the other 



hand, for a different model [22], it was found that 77 = 0.26 ±0.03 and v = 0.78 ±0.03. Both 



simulations are consistent with v e in eq. (I90p for a bosonic spinon and t 



m 



However, 



they are also consistent with a fermionic spinon with t = m. Our formulas give rj « 0.67 
and rj ~ 0.42 for fermions and bosons, respectively. 



VIII. CONCLUSIONS 

To summarize, we proposed that spin-^ particles can be described by a symplectic fermion 
quantum field theory as an alternative to the Dirac theory if one demands only rotational 
invariance rather than the full Lorentz invariance. The resulting lagrangian has a form 
resembling that of a scalar field but with the "wrong" statistics. A hidden Sp(2N) symmetry 
allows identification of spin via an SO (3) subgroup for any N, so that the statistics of the 
field is consistent with the spin-statistics connection for spin-| particles. The hamiltonian 
is pseudo-hermitian, and this is sufficient to guarantee a unitary time evolution. The usual 
spin-statistics theorem for this kind of field theory is circumvented, because the proof of the 
latter does not allow for a pseudo-hermitian hamiltonian. 

We have analyzed the renormalization-group properties and critical exponents of the 
symplectic fermion model up to two-loop order. The anomalous dimensions and /3-function 
of the Sp(2N) model are related to those of the O(M) Wilson-Fisher model by setting 
M = —2N. This correspondence between O(M) and Sp(2N) models does not hold for all 
physical properties however. In addition to the usual exponents, we have also computed 
exponents for fields that are bilinear in the fundamental spinon fields. 

The potentially most interesting possible applications of our theory are to quantum crit- 
ical spinons in d — 2 spacial dimensions. We have computed the critical exponents for 
"magnetic" order parameters that are quadratic in the spinon fields, as in models of decon- 
fined quantum criticality. Comparison with existing numerical simulations is to some extent 
favorable, but not yet conclusive. 
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